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ABSTRACT 

We present a refined orbital solution for the components A, B, and C of the nearby late-M type multiple system LHS 1070. By 
combining astrometric datapoints from NACO/VLT, CIAO/SUBARU, and PUEO/CFHT, as well as a radial velocity measurement 
from the newly commissioned near infrared high-resolution spectrograph CRIRES/VLT, we achieve a very precise orbital solution 
for the B and C components and a first realistic constraint on the much longer orbit of the A-BC system. Both orbits appear to be 
co-planar. Masses for the B and C components calculated from the new orbital solution (M B+ c = 0.157 ± 0.009 Mq) are in excellent 
agreement with theoretical models, but do not match empirical mass-luminosity tracks. The preliminary orbit of the A-BC system 
reveals no mass excess for the A component, giving no indication for a previously proposed fourth (D) component in LHS 1070. 
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1. Introduction 

LHS 1070 (GJ 2005, LP 881-64) is a low -mass triple system at 
a distance of 7.72 + 0.15 pc from the Sun dCosta et al.U2005l) . It 
consists of two components with masses close to the substellar 
limit (LHS 1070 B and C) with spectral types of M8.5 to M9, 
orbited by a slightly more massive M5.5 type component, LHS 
1070 A. 

This system is important when defining the empirical mass- 
luminosity relationship at the low-mass end of the main se- 
quence. Few reliable and accurate mass determinations, such as 
dynamical masses from binary and multiple systems, exist be- 
low 0.1 Mq, which makes LHS 1070 an important anchor point, 
especially given its age of several Gigayears and its proximity to 
the Sun. 

iLeinert et al.l d200ll) have presented the first orbital solution 
for the LHS 1070 BC system based on speckle measurements 
and early adaptive optics images. Deriving a combined mass of 
0.138 M© for LHS 1070 BC, they put the system in the con- 
text of the empiric al and theoretical m ass-luminosity relation- 
ships presented in ILeinert et al.l (|2000|) . Given the much longer 
orbital period of the outer, more massive component LHS 1070 
A, mass estimates for LHS 1070 A and a first reliable solution 
for the outer orbit were beyond reach in 2001. Nevertheless, 
ILeinert et alj d200ll) succeeded in giving rough constraints on 



* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 60.A-9078(A) and 79.C-0106(A), 
as well as data collected at Subaru Telescope, which is operated by 
the National Astronomical Observatory of Japan. Based also on ob- 
servations obtained at the Canada-France-Hawaii Telescope (CFHT) 
which is operated by the National Research Council of Canada, the 
Institut National des Sciences de lUnivers of the Centre National de 
la Recherche Scientique of France, and the University of Hawaii. 



the outer orbit, suggesting co-planarity with the inner orbit of 
LHS 1070 BC. 

Late-type double and multiple systems are also important 
probes for the connection of rotation and magnetic activity at 
the transition from partly to fully convective stars. During the 
scien ce verification phase of CRIRES at the VLT (IKaufl et all 
2006), we obtained spatially resolved measurements of the ro- 
tational velocity an d magnetic activity f or all three components 
of LHS 1070 (see iReiners et all 120071 for details). This mea- 
surement is essentially free of the uncertainties induced from 
the otherwise unknown inclination angle when comparing rota- 
tion and magnetic activity for objects of different spectral types, 
since the rotation axes of the components in this triple system 
are most likely spatially well-aligned. The CRIRES/VLT data- 
point also adds important information on the orbit of the triple 
system, removing the ambiguity in the true spatial orientation. 

In this paper we present our results from an improved or- 
bital fit of LHS 1070 BC and first constraints on the orbit of 
LHS 1070 A around the barycentre of the BC subsystem. Our 
observations, as well as extensive archival data, are presented in 
Section 2. The orbital solutions and masses derived from these 
data are discussed in Section 3. A summary of the most impor- 
tant results is given in Section 4. 

2. Observations and data reduction 

In addition to th e astrometric datapoints published in 
ILeinert et al.l d2001l) . we used archived images from the adap- 
tive optics (AO) cameras PUEO/CFHT and NACO/VLT cover- 
ing a timespan of nearly 6 years, see Table Q] for details. We 
also used the first s patially resolved, h igh-resolution spectrum 
of LHS 1070 ABC dReiners et all 120071) and determined radial 
velocities (RV) for the relative motions of LHS 1070 B and C, as 
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Table 1. Observation log and obtained relative positions for LHS 1070 C (Xc, Yc) with respect to component B and for the bary centre 
of LHS 1070 BC (X BC , Y BC ) with respect to component A. 



Date 


Instrument 


Pixel scale 


Filter 


Prog.ID. 


P.I. 




Y c 




Y BC 






(mas/pixel) 








(mas) 


(mas) 


(mas) 


(mas) 








Observations retrieved from archives 










2000-08-16 


PUEO/CFHT 


34.9 


K' 


00AD99 


Forveille 


25.0 


229.1 


1520 


206 


2000-12-12 


PUEO/CFHT 


34.9 


K 


00BH12 


Roddier 


-24.6 


218.5 


1526 


224 


2002-07-24 


PUEO/CFHT 


34.9 


K' 


02AF43 


Beuzit 


-265.2 


126.7 


1541 


315 


2002-09-11 


PUEO/CFHT 


34.9 


H, K' 


02BF27 


Beuzit 


-281.8 


116.9 


1535 


326 


2002-10-21 


NACO/VLT 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-295.8 


109.3 


1537 


334 


2002-11-15 


NACO/VLT 


13.24 


NB2.17 


70.C-0476(A) 


Bailer- Jones 


-303.3 


103.9 


1540 


337 


2002-12-16 


NACO/VLT 


13.24 


J, Ks 


70.C-0704(A) 


Leinert 


-314.3 


97.3 


1538 


342 


2003-01-25 


NACO/VLT 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-327.1 


89.2 


1536 


344 


2003-05-24 


NACO/VLT 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-362.4 


59.0 


1534 


369 


2003-06-11 


NACO/VLT 


13.24 


Ks 


71.C-0262(A) 


Bailer- Jones 


-365.9 


55.6 


1533 


370 


2003-06-17 


NACO/VLT 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-367.2 


54.4 


1534 


372 


2003-06-27 


NACO/VT T 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-369.8 


52.5 


1534 


375 




NAfOM T 

INrtv-VJ/ V LI 


13.24 


Ks 


71.C-0262(B) 


Bailer- Jones 


-^R(S 

JOU.U 


^7 ^ 




^78 


2003-09-14 


NACO/VLT 


13.24 


J, Ks, NB2.17 


70.C-0476(A) 


Bailer- Jones 


-388.3 


35.1 


1533 


381 


2003-12-10 


NACO/VLT 


13.24 


NB1.64 


72.C-0570(A) 


Beuzit 


-405.8 


14.8 


1532 


390 


2003-12-12 


NACO/VLT 


13.24 


J, H, Ks 


72.C-0022(A) 


Leinert 


-406.3 


14.8 


1531 


391 


2004-12-11 


NACO/VLT 


13.24 


Ks 


74.C-0637(A) 


Kohler 


-443.6 


-73.3 


1511 


446 


2006-10-30 


NACO/VLT 


13.24 


NB2.17 


78.C-0386(A) 


Ratzka 


-347.2 


-197.1 


1458 


524 








Own observations 












2006-10-09 


CRIRES/VLT 






60.A-9078(A) 


Reiners 










2007-07-16 


NACO/VLT 


13.24 


NB 


79.C-0106(A) 


Neuhauser 


-260.7 


-223.7 


1426 


563 


2007-07-25 


CIAO/SUBARU 


10.82 


NB 


o07111 


Neuhauser 


-258.2 


-224.4 


1431 


567 



" X c and Y c are accurate to +1.5 mas 
h X BC and Y BC are accurate to ±6 mas 



well as for the motion of LHS 1070 A around the barycentre of 
LHS 1070 BC. Finally, in July 2007 we observed LHS 1070 in 
two imaging campaigns with NACO/VLT and CIAO/SUBARU 

All AO imaging data were reduced by combining individ- 
ual exposures (shift and add) and determining the positions and 
fluxes of each objec t (including their errors) with IDL/starfinder 
dDiolaiti et all 120001) . a program based on empirical PSF fitting. 
In cases where a high number of frames were recorded, the er- 
rors were calculated from statistics over the distribution of the 
single positions obtained from the individual frames instead of 
the formal error computations on a single co-added frame. 

The pixel scale and detector orientation of PUEO/CFHT 
were determined from reference images of the HIPPARCOS 
double star HIP482 (p=7.687", 0=148.29°) obtained during 
each of the respective runs. The pixel scale was stable at 
34.9+0.1 mas/pixel, but the detector orientation was changing 
between the runs by up to +2.5°, and was stable during each 
of the runs on the 0.05° level. For all NACO/VLT observations 
taken with the S13 camera, a fixed pixel scale of 13.24+0.05 
mas/pixel, and a detector orientation of 8 — 0.0 + 0.1° was 
adopted. Our exper ience from othe r long- term monitoring pro- 
grammes (see i.e. iNeuhauser et all 120071) show that both, the 
pixel scale and the detector orientation of NACO/VLT are very 
stable, and the accuracy of their absolute value is limited by the 
astrometric errors of available HIPPARCOS double stars rather 
than instrumental instabilities. Finally, our CIAO/SUBARU im- 
ages were calibrated by comparing an HST/WFPC image of the 
core region of the globular cluster Ml 5 with images of the very 
same region in M15 taken during our run with CIAO/SUBARU. 



We note that for the NACO/VLT observations the dominant 
error source is the absolute pixel scale and detector orientation 
rather than the centroiding error on the targets. The latter can 
easily be below 1 mas, especially when enough frames with good 
S/N were obtained. However, the typical accuracy of NACOs 
S13 pixel scale is about ±0.05 mas/pix, which in turn limits the 
accuracy of the determined positions to about +1.5 mas for the 
typical separation of LHS 1070 BC (~ 0.5") and +6 mas for the 
much larger separation of LHS 1070 BC around component A. 
Datapoints being closer in time than a few days were combined 
to a mean epoch since the expected motions are far smaller than 
the derived errors. 

The data reduction and calibrat i on of our CRIRES/VLT 
spectrum is outlined in lReiners et al.l d2007). Since here we de- 
rive relative radial velocities, our result critically depends on a 
good dispersion solution. The covered spectral range falls into 
a clear atmospheric window. CRIRES is equipped with a Th-Ar 
hollow cathode lamp for wavelength calibration dKerber et all 
120071) . More than 2000 lines have been established as wave- 
length calibration standards by laboratory measurements in a 
collaboration between the US National Institute of Standards 
and Technology (NIST) and ESO (Kerber et al. in preparation). 
These calibration reference data are implemented in the CRIRES 
data reduction package. We made use of this linelist to calibrate 
our data. 

For the respective relative radial velocity between each of 
the components of LHS 1070, the alignment of the components 
in the slit is crucial. At the time of the observation, all three 
components are nearly, but not perfectly, aligned along the slit. 
The misalignment (as a deviation from a perfect linear configu- 
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ration) was ~9°. From geometrical considerations (and by taking 
the individual S/N ratios into account), we adopt +1.5 km/s for 
the uncertainty in the relative RV between LHS 1070 B and C, 
and for the relative RV between LHS 1070 A and the barycen- 
tre of LHS 1070 BC. We note that the most important constraint 
for the orbital solution is given by the sign of the radial velocity 
rather than by its exact value. 

3. Analysis and results 

3.1. Astromethc solution and orbital fitting: LHS 1070 BC 

To fit the relative orbit of LHS 1070 B and C around their com- 
mon bary centre, we chose a similar approach to the one pre- 
sented in iLeinert et al.l d200ll) . We simultaneously s olved the 
Kepler ian equation using an IDL routine based on iHilditchl 
(2001) for all astrometric and spectroscopic datapoints. A x 2 
minimisation was performed to find the optimal se t of pa ram- 
eters, and a Jacknife approach dShao & Dongsheng, 1996) was 
chosen to determine the errors in the orbital parameters (see 
Table [2] for the parameters and Figure Q] for the orbital fit). The 
orbital coverage is now 253°. The reduced ;f 2 1.27 is rea- 
sonable but points towards a slight un derestimation of the real 
uncertainties in the original dataset of ILeinert et al.l d200ll) and 
in the NACO/VLT data. 

We note that the refined orbita l parameters of the B C system 
agree with the ones published in ILeinert et alj ([2001 ), but only 
the time of periastron (7q) differs slightly, being 1991-02-18 in 
our calculation. The error on the orbital period and semimajor 
axis is now at the sub-percent level, allowing a very precise mass 
determination for a fixed-distance estimate to the system. 

Table 2. Orbital and system parameters for LHS 1070 B and C. 



Parameter 



Value 



Error 



«(") 


0.4619 


+0.0007 


e 


0.034 


±0.002 


in 


62.59 


±0.12 


no 


14.50 


±0.12 




240.33 


±0.24 


P (days) 


6214.74 


±0.42 


T (MJD) 


48305.6 


±1.8 


Ktrig (") 


129.47^ 


±2.48 


d (pc) 


7.72 ~ 


±0.15 


a (AH) 


3.57 


±0.07 


M B + M C {M Q ) 


0.157 


±0.009^ 



" from lCosta etafl d2005h 

b including error in distance, at fixed distance the error is 0.0007 Mq 
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Fig. 1. Relative orbit of LHS 1070 C ar ound component B . 
Upper panel: astrometric observations from lLeinert et al.l (1200 ll) 
(black crosses) including error bars, PUEO/CFHT archival dat- 
apoints (blue asterisks), NACO/VLT archival datapoint (red di- 
amonds), along with a new NACO/VLT datapoint (green box) 
and a new CIAO/SUBARU datapoint (blue hexagon). Errors and 
deriv ations (O-C) for all datapoints except the original points 
from Lei nert et al.l d2001l) are smaller than their respective sym- 
bols. The time of the radial velocity measurement (2006-10-09) 
is marked with a circle, the periastron with a triangle. The line of 
nodes is shown. Lower panel: Relative radial velocity between 
component C and B. The CRIRES/VLT datapoint is shown with 
its lcr error bar. 



3.2. Astrometric solution and orbital fitting: LHS 1070 A 

The datapoints for the orbit of the BC system around LHS 1070 
A presented in lLeinert et al.l ([2001 ) did not show any strong cur- 
vature, hence, the authors used a general x 1 fit with combina- 
tions of period and eccentricity to map the parameter space. They 
used a stability argument for hierarchical triple systems to fur- 
ther constrain the period and inclination angle. The new data- 
points reveal a strong curvature, the ascending node being prob- 
ably close in time to our CRIRES/VLT radial velocity datapoint. 
Thus, the new datapoints give a stronger constraint on the outer 



orbit, event hough the covered arc of the orbit is still only about 
33°. However, the orbital solution is still not unique. We per- 
formed the^ 2 minimisation for fixed pairs of eccentricity and or- 
bital period, spanning a grid of e — 0...0.45 and P = 40... 120 yr. 
Three families of solutions appeared to fit the data: first a small 
region where the eccentricity is zero and the orbital period is ~ 
82 years, and second, two narrow regions for rising eccentricity 
with either shorter or longer periods with the longitude of peri- 
astron (o>) turning by 180° from ~ 150° to ~ 327° between these 
two solutions (see Fig. [3). 



4 



A. Seifahrt et al.: Improved orbital solution and masses for LHS 1070 



1.5 



L 2007-07-25 



u 

w 
Q 

>< 



1.0 r 



0.5 r 



0.0 r 



-0.5 r 



-1.0 r 



1.5 




1993-07-29 



0.5 0.0 -0.5 
Y: RA (arcsec) 



-1.0 



-1.5 



1.4 



u 
Q 

X 



1.2 



= = 








L a 






























X. j 










0.6 


0.4 


0.2 


0.0 -0.2 



Y: RA (arcsec) 



% 

i 
< 



4 ir- 

2 =- 

^ 
-2 T 
-4 =- 





5.0-10 3 1.0-10 4 1.5-10 4 2.0-10 4 2.5-10 4 3.0-K 

MJD - 50537.3 



Fig. 2. Relative orbit of the barycentre of LHS 1070 B and C 
around component A. Upper panel: astrometric observations. 
Same notation and symbols as in FigureQ] Central panel: Cutout 
of the upper panel for better visibility. Lowest panel: Relative ra- 
dial velocity between the barycentre of LHS 1070 B and C and 
component A. The CRIRES/VLT datapoint is shown with its lcr 
error bar. 



We have to use a st ability criterion to f urther constrain the 
parameter space, as in Lei nert et al.l ( 20011). Here we use the 
critical distance parameter ^/^l fromlHarringto n] (| 19751). where 
2 = a %{\ —e%), with updated values from lDonnison & Mikulskia 
1995b to infer the stability limit of the LHS 1070 system. Given 
a prograde orbital motion, the mass estimates of the three bodies 
and solutions for the semimajor axis of the BC system (a\), as 
well as the values for the eccentricity (ej) and semimajor axis 
(02) of the wide orbit from our grid solution, we can plot the 
line of stabilitjQ {qila\ > 3.3) in the P-e plane of Fig. [3] Only 
solutions below and to the right of this line are long-term sta- 
ble, excluding most of the regions with orbital periods shorter 
than 80 years and eccentricities greater than zero. This is similar 



1 Stability is here defined as no significant change in the orbital pa- 
rameters over a significant number of orbital revolutions. 
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Fig. 3. The most likely values of eccentricity and orbital period 
for the outer orbit (barycentre of LHS 1070 BC around compo- 
nent A). The contours show the 1,2, and 3<x contours around the 
X 2 minimum. The left branch represents a family of solutions 
where a> 150°, while the right branch is for a family of solu- 
tions where u> 327°. The dash-dotted line marks the stability 
threshold of the system (see text for details). Only solutions be- 
low this line are long-term stable. 



to what lLeinert et al.l d2001l) propose, based on slightly different 
assumptions. 

Having narrowed down the phase-space of possible solutions 
for the outer orbit considerably, we evaluate the histograms of 
the most important orbital parameters for solutions within 2<x 
around min(^ 2 ) and within the stability region (see Fig. |4). We 
find that the distribution of the inclination angle for the outer 
orbit is within 2° of the value found for the BC orbit (~ 62.6°) 
and D. is nearly identical for both orbits (~ 14°). Hence, the inner 
and outer orbits appear to be co-planar. In Fig. [2] we show the 
datapoints and the most likely orbital solution with P ^ 82 yr 
ande=0.01. 



3.3. Masses 

With the new orbit and refined distance dCosta et all 120051) of 
the LHS 1070 system, we can constrain the combined dynamical 
mass of the B and C components to M B +M C = 0.157+0.009MQ. 
The dominating error source for the dynamical mass of LHS 
1070 BC is still the uncertainty of +0.15 pc in the distance of 
the system. At a fixed distance of 7.72 pc, the error in mass 
is as small as 0.0007 Mq. Our value for the combined mass 
of LHS 1070 BC i s higher by 9% than the one obtained by 
iLeinert et al.l (12 001 ). mainly because the new distance estimate 
by ICosta et aU (|2005[) is highe r by nearly 3% than the value 
adopted by ILeinert et al.l (12001 1) . This improvement moves the 
da tapoints for LHS 10 70 B and C closer to the 5 Gyr isochrone 
in ILeinert et al.l d200ll Fig. 5); see also Fig. [5] here. The com- 
bined dynamical mass of LHS 1070 B and C is thus in very good 
agreement with the model predictions for the com bined mass 
of LH S 1070 B a nd C of 0.159 . 0.1 61 M from ILeinert et al.l 
(120001) (based on lBaraffe et~aLl Il998l) . Hence, both components 
are slightly above the hydrogen-burning mass limit. 

We note that the mass predictions from empirical mass- 
luminosity re lations differ s from our value by a consider- 
able amount. ILeinert et al .] (120001) give a combined mass for 
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il (deg) Mass of LHS 1070 A (MJ 

Fig. 4. Histograms of the distribution of orbital parameters from 
our grid of possible solutions for the outer orbit within 2<x 
around min(^ 2 ) and within the stability region inferred for the 
outer orbit. The mass of LHS 1070 A was derived from the com- 
bined system mass and the mass of LHS 1070 BC (see Table 1). 
The spread in mass does not include errors from the distance of 
the system. The dashed-dotted lines indicate the respective val- 
ues from the inner orbit (LHS 1070 BC) for comparison. 




0.5 0.4 0.3 0.2 0.1 
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Fig. 5. Mass-luminosity r elationship at the end of the main se- 
quence. Data taken f r om | Delfosse et alj d2000l) . except for GJ 
1245 C dHenrv et all 1 19991) a nd LHS 1070 ABC (this paper, 
My from iLeinert et all l2000h . The solid l ine shows t h e em- 
pirical mass-luminosity relationship from iHenrv et alj (|l999, 
Equ. 7) and the da shed line the theoretical relationship from 
Baraff e etal.l(ll998l) . 



tracks of 0.097...0.113 M© (ILeinert et al.L 120001; iBaraffe et all 
120031) . 

In Fig. [5] we show the mass- lumi nosity relation with data- 
points from lDelfosse et al.1 (l2000h and lHenrv etaTJ d!999h . LHS 
1070 ABC is shown in this graph with absolute V magnitudes 
obtained from lLeinert et al.l y000) and our new mass estimates. 
Note that the error margin for the mass of LHS 1070 B and C 
are shown twice as large to accommodate the uncertainty in the 
mass ratio of both components. As becomes apparent from this 
plot, the empirical mass-luminosity relation is too flat, underpre- 
dicting masses belo w 0.1 M©, while the model prediction from 
Baraffe et al. ( 1998) is much steeper and fits the available data- 
points. 



LHS1070 B and C of =* 0.138M o from optical, and =* 
0.1 15...0.130Mq from the ne a r infr ared mass-luminosity re- 
lations based on IHenrv et all (1 19991) and IHenrv & McCarthy 
(1993). Such predictions underestimate the true mass by 12 to 
26%. 

The first constraint on the orbit of LHS 1070 A around com- 
ponents BC allows the calculation of a combined mass for 
the whole system. The error budget for the combined mass of 
LHS 1070 ABC is still dominated by the uncertainties in the pa- 
rameters of the outer orbit (semimajor axis and orbital period). 
However, orbital period and semimajor axis are strongly corre- 
lated, and their respective error contributions do not simply add 
up. The combined dynamical mass of the whole system of LHS 
1070 ABC is M LHS i070 = 0.272 ± 0.017M©. 

Subtracting the well-constrained dynamical mass of 
LHS 1070 BC from the total dynamical mass of the system 
leaves M A = 0.115 + 0.010M© for the A component. The 
precision is higher than for the total mass estimate of LHS 
1070 ABC since the error contribution from the uncertainty in 
distance nearly cancels out. The resulting dynamical mass for 
LHS 1070 A is only slightly higher than the predictions, both 
from theoretical models of 0.098. ..0.109 M© and the empirical 



3.4. Is there a fourth component in the system? 

IHenrv et"aTI (1 1999b report the discovery of a fourth component 
in the LHS 1070 system from HST FGS measurements. The au- 
thors give a magnitude difference of AV 2.5 to LHS 1070 A, 
a mass of 0.07 M© for LHS 1070 D, and a separ ation of ~ 50 
mas fr om LHS 1070 A at their HST measurements. lLeinert et all 
(2000) find no clear presence of such a rather massive, but 
still much cooler and unresolved component in the spectrum of 
LHS 1070 A. 

We checked all images, especially the ones observed with 
NACO/VLT under good atmospheric conditions (hence with a 
good strehl ratio) for any additional object close to LHS 1070 A. 
We did not find any object at distances equal to or greater than 
the diffraction limit of the VLT at J, H, and Ks band (32 to 
56 mas, respectively). We did not find any indication of an elon- 
gation or any difference in the PSF shape of LHS 1070 A, B, 
and C. This rules out the existence of a component with more 
than 0.060 M© at distances over 32 mas of AJ < 4.5mag and at 
distan ces over 56 mas of AKs < 5.0mag, based on Baraffe et al. 
(2003) and our detection limits. Given the extended time cover- 
age and tight sampling of the NACO/VLT datapoints, it is rather 
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unlikely that we have always missed the object because the dis- 
tan ce of both componen ts was too small to be resolved. 

iLeinert et al.l d200ll) argue that the reflex motion of an ob- 
ject of intermediate mass between LHS 1070 A and D should 
induce a reflex motion of the order of 10 mas, which was not de- 
tectable with their dataset. Analysing the residuals of the later 
NACO/VLT measurements, we indeed find a scatter of of 1- 
3 mas, exhibiting a non-Gaussian distribution. However, given 
the large separation of the observations in comparison to the pre- 
sumed orbital period of the tentative fourth component, we are 
unable to find an orbital solution for these points. We can only 
argue that the detec ted derivations are considerably smaller than 
the prediction of ILeinert et al.l (1200 ll) . implying a much lower 
mass than previously anticipated. 

Our mass estimate of LHS 1070 A is in good agreement with 
the empirical and theoretical predictions. We do not find any in- 
dication of a fourth component in the triple system LHS 1070 
with mass down to ~10 Mj up . 

4. Summary 

Using a combination of published and archival data and adding 
new observations from high spatial and spectral resolution imag- 
ing and spectroscopy, we have considerably improved the orbital 
solution for the low-mass triple system LHS 1070. Most notice- 
ably we found that: 

- The dynamical masses for LHS 1070 B and C, both close to 
the hydrogen burning limit, are in good agreement w i th the- 
oretical predictions from the models of Bar affe et al.l dl998l) 
but higher by up to 25% than the empirical near infrare d 
mass-luminosity relationship of lHenrv & McCarthy! (ll 993b . 
Since only one datapoint (GJ 1245 C) is constraining the 
mass-luminosity function below My = 18 (see Fig.|5]l, the 
addition of an accurate value for LHS 1070 BC marks an im- 
portant anchorpoint for the mass-luminosity relation at the 
end of the main sequence. 

- The orbits of the LHS 1070 BC system and the outer orbit of 
LHS 1070 A around the barycentre of LHS 1070 BC appear 
co-planar: their inclination angle agrees to within 2° and the 
respective longitudes of the ascending node (£2) of both or- 
bits are identical within their uncertainties. 

- The mass estimate of LHS 1070 A matches the predictions 
from theory and from the empirical tracks. Given its error 
margin and the fact that no companion could be identified in 
the well-resolved VLT/NACO images, we find no indication 
for a fourth component (LHS 1070 D) in close orbit around 
LHS 1070 A. 
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